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Abstract 

The order of the vortex state in Lai.9Sro.iCu04 is probed using muon spin rotation and smah- 
angle neutron scattering. A transition from a Bragg glass to a vortex glass is observed, where 
the latter is composed of disordered vortex lines. In the vicinity of the transition the microscopic 
behavior reflects a delicate interplay of thermally-induced and pinning-induced disorder. 

PACS numbers: PACS numbers: 74.25.Qt,74.72.Dn, 76.75.+i, 61.12.Ex 
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It has been acknowledged for many years [l| that the long range translational order of 
the Abrikosov vortex crystal in type II superconductors should be disrupted by arbitrarily 
small amounts of disorder. The current theoretical perspective is that weak disorder results 
in a 'Bragg glass' (BG), which retains the topological order of the flux line lattice (FLL) 
but yields broadened diffraction peaks [2]. Strong disorder should result in the topologically 
disordered 'vortex glass' (VG) [2]. The BG concept accounts for the observation of an 
apparently ordered FLL in many systems (e.g. 0, 4]) and further evidence for this phase in 
an isotropic superconductor has recently been obtained using small angle neutron scattering 
(SANS) ^. At high magnetic fields and low temperatures, where the effects of disorder are 
stronger, the BG should be unstable to the formation of the VG. 

In contrast to the enormous wealth of experimental work on vortex matter derived from 
macroscopic measurements , |9|, there have been relatively few microscopic studies. 

Of these techniques SANS and muon spin rotation (uSR) measurements provide unique 

n n fi n 

information concerning the local and long range order of the vortex system (e.g. [3|, U, I5l.ll0|). 
Here we report muon-spin rotation measurements on underdoped Lai.9Sro.iCu04 (LSCO) 
which provide the first unambiguous experimental evidence for the transition from an ordered 
phase (BG) to a VG in a system of well-coupled vortex lines. 

The importance of LSCO as a system in which to study vortex matter, and in particu- 
lar the BG to VG transition, can be understood in terms of two key parameters, namely 
the superconducting penetration depth Xab and the anisotropy parameter 7. The former 
determines the ability of the supercurrents in the copper-oxide sheets (a6-planes) to screen 
externally applied magnetic fields, while the latter quantifies the uniaxial anisotropy of the 
penetration depth. In addition, the Josephson length Aj = 7s (where s is the spacing of 
the copper oxide planes), determines the effectiveness of currents, tunnelling between the 



conduction planes, in maintaining the stiffness of a vortex line . While higher values 
of Xab, 7 increase the susceptibility of the vortices to disordering, it is the ratio Xj/Xab 
which determines the dimensionality of the vortex system. For example, in optimally doped 
YBa2Cu307_5 (YBCO) Xj/Xab ^ 1, so the vortices resemble the rigid rods of flux found in 
conventional isotropic superconductors. In contrast, for optimally- doped Bi2Sr2CaCu208+5 
(BSCCO) Xj/Xab ~ 1; consequently the vortex system is quasi-two dimensional. As we 
demonstrate here, the novelty of underdoped LSCO is that the material parameters of 
{Xab ~ 3000 A and Aj ~ 250 A ) give rise to a system of fairly rigid vortex lines which are 
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nonetheless highly susceptible to transverse fluctuations. 

Our fiSK measurements were performed using the GPS spectrometer on the vrMS beam- 
line at the Paul Scherrer Institut, Switzerland. The fleld was applied parallel to the c-axis of 
the crystal, and the muon spins were initially polarised nearly perpendicular to the applied 
fleld. The muons implant in the sample at random positions with respect to the vortices, 
and each precesses at a rate u = '-)^B{r), determined by the local fleld B{r) and the gyro- 
magnetic ratio of the muon. The muon spin precession is monitored via the positrons which 
are emitted preferentially along the muon-spin direction during the decay of the muon lo| . 
The precession signal, averaged over ~ 10^ muons, reveals the probability distribution p{B) 
of the spatial variation B{r). The sample was a single crystal of length 20 mm and diam- 
eter 5 mm, with the c-axis perpendicular to the axis of the cylinder. The superconducting 
transition temperature of the sample was Tc = 29 K with a width (10% - 90%) AT = 1.5 K, 
indicating the high quality of the sample. All yuSR measurements were made after cooling 
in an applied fleld in order to obtain a uniform flux density across the sample. 

In a /iSR experiment the presence of a vortex line lattice is revealed by a p{B) which has 
a distinct signature, most notably a peak followed by a tail extending towards flelds higher 
than the mean, reflecting the small volume of the core region of the vortices. The peak 
at the mode of the distribution represents the high probability density of the fleld value 
at the centre of a line connecting two vortex cores [10]. In Fig. [1^ are shown probability 
distributions measured in LSCO at low fleld, where those taken at the lowest temperatures 
correspond most closely to the hneshape expected from an ideal lattice. One of the calculated 
lineshapes given in Fig. [lb is also derived from an ideal vortex lattice, where for clarity the 
rounding effects of instrumental resolution and nuclear dipolar broadening, visible in the 
experimental data of Fig. [1^, have been omitted. Departures from this ideal form may be 
used to deduce changes in the morphology of the vortex state [lo| . The muon probe measures 
the local fleld, which is determined by contributions from vortices only within ~ \ab of the 
muon; hence, the general form of p{B) is remarkably robust to changes in long range order 
12l |. In contrast, the technique can be highly sensitive to changes in the local environment 



lOj , which allows us to observe the 



due to thermally-induced or pinning-induced distortions 
transition from the ordered state to the VG phase in LSCO. 

The main result of this paper can be appreciated immediately by comparing the nor- 
malised experimental lineshapes of Fig. [TK and [lb, measured at 80 Oe and 6 kOe respec- 



4 



tively. Referring first to the lowest temperature lineshape at eacli field, it is clear that they 
are very different. While the 80 Oe distribution has the characteristics of an ideal vortex 
lattice described above, the 6 kOe signal is highly symmetric, indicating a strong departure 
from an ideal vortex lattice or Bragg glass arrangement. For comparison, in Fig. [T]: we 
include the p{B) derived from Monte Carlo simulations of a perfect triangular Abrikosov 
lattice and also for a vortex glass structure having short-ranged translational correlations of 
the order of about six inter- vortex spacings. The simulations are based on one possible mech- 
anism for the loss of long range order recently discussed by one of us 13|. In this scenario, 
the system undergoes a transition from a BG phase to a multi-domain glass comprising a 
size distribution of domains within which the vortices are locally ordered. Increasing the 
field leads to a rapid fall in the average domain size just above the BG phase. Note that the 
distribution function obtained for the multi-domain structure is far more symmetric as well 
as broader than is the case for the perfect lattice, as is found experimentally in the data at 
6 kOe (see Fig. [Do). 

Further conclusive evidence for the transition from a BG to a VG phase comes from 
examining the width of the lineshapes a, given by the square root of the second moment of 
p{B). This quantity is plotted in Fig. [2] as a function of applied field, for a range of temper- 
atures. We first draw attention to the data at the lowest temperatures, where it can be seen 
that the signal measured at 6 kOe is considerably broader than that at 80 Oe. It has been 
shown theoretically jl^ . fisj ] that such a broadening of the signal from the vortex lattice can 
only arise from static disorder in the positions of vortex lines within a plane perpendicular 
to the field 16| . This arises due to positional disorder of the vortices, which gives rise to re- 



gions with field values both higher and lower than in the well-ordered lattice [1^, [ij, Il5|, [ITJ ■ 



Conversely, Brandt showed that in a system composed of two-dimensional 'pancake-vortex' 



strings U, short-wavelength transverse fluetuafons along the field d.rection would always 
lead to a narrowing of the field distribution p{B) p^, and indeed this was the situation 



found experimentally in the very anisotropic (quasi-two-dimensional) BSCCO [181]. Thus 
it is clear that in LSCO we are dealing with a highly disordered vortex line arrangement, 
which we identify with the VG phase. 

We now turn our attention to the temperature dependence of the hneshapes of Fig. [1^, 
where it can be seen increasing temperature truncates the high-field tails. This reflects the 
increasing amplitude of thermally- induced fluctuations (u^(T))^/^ of the vortex positions. 
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on a timescale faster than the characteristic muon samphng time p^. The muons thus 
experience a time-averaged field distribution in which the high-field values arising from close 
to the vortex cores are smeared out. Compared to those at low temperature, the lineshapes 
are narrowed. Such narrowing arises both from the increase in A(T) and (n^(T))^/^. This 
may be described by [lol, Q : 

-r2(„2>/2 



where the f are reciprocal lattice vectors, multiples of n/d, and d ~ y^o/B is the plane 
spacing of the vortex lattice. At a given temperature, the width a decreases as the field is 
increased, mainly due to an increase in the ratio {u'^y^'^/d. This is in contrast to the case of 
an ideal static vortex lattice, where a is independent of field in the range Hc^ < H <^ Hc^- 
In Fig.[2]one observes a reduction in width with increasing field up to fields of around 1 kOe. 
While the observed changes of the lineshapes show that this must at least partly be due to 
these dynamical effects, simulations indicate that the magnitude of this narrowing cannot 
be entirely attributed to this mechanism. The broadening towards lower fields must also 
include a contribution from an increase in the static disorder due to the reduction of elastic 
interactions between vortices in this very dilute vortex state. Indeed, it is only at fields 
above a few hundred Oersted that the vortex separation becomes less than the penetration 
depth in this system, whereby the shear modulus for the lattice increases significantly 
Thus below the BG-VG transition one would also expect the linewidth in this system to 
increase with decreasing field when vortices cease to overlap significantly. A minimum in 
(j{B) thus occurs due to the competition of these broadening mechanisms in the vicinity of 
the transition to the VG phase. 

We have also recently used SANS to perform a preliminary experiment to image the 
vortex arrangement in the same LSCO sample, using instrument D22 at the Institut Laue- 
Langevin, Grenoble, France. At an applied field of 150 Oe an hexagonal diffraction pattern 
could be measured (Fig. [3^), thus confirming the quasi-long-range order of the vortex lattice 
at low fields. At a temperature of 5 K the intensity was found to drop very rapidly as a 
function of applied field. Furthermore, the diffraction spots show an appreciable broadening 
in q in the vicinity of Bcr ~ 800 G 21|, which, in accord with the /iSR data, implies an 
increase of static disorder in the system around this field. While in /iSR experiments a 
field-cooled (FC) experiment is necessary to interpret the lineshapes, in other systems this 



is known to not always represent the most ordered state 22]. In the SANS experiment 
we also attempt to improve the lattice perfection in the Bragg-glass region by cooling in a 



modulated applied field, as in ref. [22]. The lattice perfection was always found to be reduced 
by these approaches, indicating that the FC state for both SANS and fiSK represents the 
most ordered vortex lattice state at low fields. 

In the literature, evidence for the existence of the BG-VG transition is frequently taken 
from signatures in the bulk magnetisation, since the transition to the glassy state is accom- 
panied by an increased effectiveness of point pinning. In Fig. [3]d is a section of an M{H) 
loop at 16.5 K, where a 'peak effect' is evident as an anomalous increase in the magnitude of 
the magnetisation beginning at around 400 Oe. The temperature dependence of the position 
of the onset of this peak i?o„(T) is plotted in the phase diagram of Fig. HI Similar features 
have previously been associated with the BG-VG transition e.g. [23i^. This peak disappears 
at around 25 K, as indicated by the line Bon(T) of Fig. HI the same temperature at which 
Bcr{T) intersects the irreversibility line Birr{T) (Fig. Hj). 

The occurrence of an equilibrium phase transition with increasing magnetic field has been 
discussed theoretically by many authors includingGiamarchi and Le Doussal It is worth 
noting in this context that the muon is sensitive only to the local magnetic field, so in general 
this transition will be manifest as a crossover of behaviour which refiects the underlying sharp 
transition. In Fig.lHthe line Bcr{T) is a plot of the minima in cr{B), where the broadening due 
to static glassy disorder begins to dominate. This feature thus provides an upper limit for the 
BG-VG transition. An upward trend of this BG-VG transition, reminiscent of that observed 
experimentally, -8^(7"), has been predicted to arise from the reduced effectiveness of pinning 
at higher temperature 2|,|2J]. The dependence of the signature observed in the magnetisation 
measurements, i?o„(T), is very differerent to Bcr{T). This is not surprising given that 
the muons measure p{B) in a field-cooled state, whereas the magnetisation measurements 
determine the macroscopic properties of a disturbed system possessing strong macroscopic 
fiux gradients, and refiect the changing dynamic response to the pinning landscape. The 
SANS intensity, by contrast, is sensitive both to thermal fiuctuations and to the loss of long 
range coherence in the vortex structure. 

In conclusion we have measured the vortex state at a microscopic level and observed a 
change with field from a BG to a VG phase, the latter possessing only short range transverse 
order. The minimum in the fiSR. linewidth with increasing field, observed at Bcr{T), arises 
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primarily from the competing influences of increasing elastic interactions combined with 
thermal fluctuations, and a rapid collapse of translational correlation lengths in the vortex 
system within the plane perpendicular to the field direction. This provides unambiguous 
evidence for a transition with increasing field from a nominally Bragg glass phase to a more 
disordered vortex glass state. This is the first such measurement on a system composed of 
well-coupled vortex lines, and paves the way for a detailed study of the evolution of short 
range order in the presence of weak pinning, a problem of universal significance j^. 
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FIG. 1: (Color) a) Temperature-dependent probability distributions of the internal field, p{B), 
due to the presence of the flux lines at applied fields of a) 80 Oe and b) 6 kOe. The curves 
are normalised at the mode of the distribution to allow easy comparison of the shapes, c) Field 
distributions derived from Monte Carlo simulations of the vortex state at 6 kOe (see text) for the 
case of: i) an ideal vortex lattice ; ii) a VG phase with a transverse correlation length of about six 
lattice spacings. The latter also includes a small contribution due to the effects of thermal narrowing 
(as observed experimentally). The curves are again normalised to the mode of the distribution. 
Note the increase in width and symmetry of the lineshape due to the increased disorder. 
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FIG. 2: (Color) a) The variation with field of the line width a at selected temperatures. The initial 
reduction with increasing field corresponds to the field dependent influence on the field distribution 
of both thermal fluctuations and increasing elastic vortex interactions. Above Bcr ~ 1 kOe, the 
width dramatically rises due to the onset of static disorder in the VG phase. Note that the high- 
field dependence has almost disappeared by a temperature of 23 K, and disappears completely at 
the irreversibility line (see Fig. H]) . 
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FIG. 3: (Color) a) A small angle neutron diffraction pattern showing the existence of an ordered 
vortex phase at low fields (150 Oe) and temperature (6.2 K), having hexagonal symmetry. The 
intensity of this pattern rapidly falls as the field is increased towards the VG phase (see text), 
b) The feature in the magnetisation at 16.5 K is associated with the crossover to the VG state, 
measured on a small piece taken from the fiSH sample. The inset shows a complete hysteresis loop. 
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FIG. 4: (Color) The magnetic phase diagram derived from the changes observed in the ;uSR field 
distributions. Bcr{T) indicates the onset of the broadening at high field, which is significant only 
below ~ 25 K, and should be considered as an upper limit for the BG-VG transition (see text). 
This uncertainty in the exact position of the transition is represented schematically by the shading 
below the line Bcr (T) . Birr (T) indicates the irreversibility line as determined by bulk measurements 
of the field-cooled/zero-field-cooled (FC/ZFC) magnetisation using a SQUID magnetometer. The 
position of the feature in the magnetisation curves Bon{T) is also plotted (see FigJSja). 
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